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Abstract 



We extend the spectral theory used for the calculation of the effective 

linear response functions of composites to the case of a poly crystalline material 

with uniaxially anisotropic microscopic symmetry. As an application, we 

y^O . combine these results with a nonlinear decoupling approximation as modified 

by Ma et al, to calculate the third-order nonlinear optical susceptibility of 

a uniaxial polycrystal, assuming that the effective dielectric function of the 

polycrystal can be calculated within the effective-medium approximation. 
OV 
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I. INTRODUCTION 

Almost twenty years ago, Bergman |I| developed the spectral approach for calculating the 
dielectric constant and other linear response functions of a two-component composite. His 
approach was to study the analytical properties of the effective dielectric constant, viewed as 
a function of the ratio of the dielectric constants of the constituents. Among other results, 
he showed that all poles of this function can be expressed as eigenvalues of a certain linear 
boundary- value problem, while the residues of those poles are given as certain integrals 
over the corresponding eigenfunctions. Bergman's theory actually describes a wide class of 
mathematically similar physical problems in which a divergence-free field appears as a linear 
response to the gradient of a potential; thus, it can be used to find the effective electrical 
and thermal conductivities, magnetic permeability and many other effective parameters 
described by mathematically equivalent equations. 

In some microgeometries, this eigenproblem can be solved by expanding the exact eigen- 
functions of the composite in terms of the individual grain eigenfunctions. This approach 
has been used to calculate the effective parameters of granular composites corresponding 
to several different ordered microgeometries [|TJ. However, the spectral representation is 
often useful even in composites where the microgeometry is not known exactly; in such 
cases, one must generally resort to various approximations in order to calculate the relevant 
spectral functions. Moreover, the spectral approach is not limited to problems involving 
linear response. For example, some of the effective macroscopic nonlinear response func- 
tions can be expressed in terms of linear response functions of the composite, and certain 
geometric factors || . Recently, this connection has been employed || , together with certain 
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approximations for linear composites, to use the spectral representation in calculating these 
nonlinear response functions. 

The present paper is directed towards the linear and nonlinear dielectric response of 
polycrystalline materials. Such materials are not generally thought of as composite media, 
but in fact they behave like composites. The reason is that even though each crystallite 
is made of the same material, it has a different spatial orientation and hence has different 
constitutive properties referred to axes fixed in the lab coordinate system. In particular, in 
this paper we extend the spectral theory to describe both the linear and nonlinear response 
of a polycrystalline material. 

By a polycrystal, we mean a material with anisotropic transport properties, such that 
the crystal symmetry axes vary in direction from point to point in space. Several previous 
workers (see, for example, JJ) have described polycrystals as composite materials. In the 
present work, we further restrict our discussion to polycrystals of uniaxial materials. In 
this case, two of the three principal components of the dielectric tensor are equal, and it 
is more straightforward to develop a spectral representation for the effective properties. 
This restriction to uniaxial materials still leave many classes of crystalline materials open to 
study. In particular, the theory should satisfactorily describe such classes of materials as the 
quasi-one-dimensional organic conductors ||, or the quasi-planar or CU2O -based high-T c 
superconductors. 

We will use this approach not only to describe the linear properties, but also to calculate 
the enhancement of the third-order nonlinear susceptibility of a polycrystal. Although this 
enhancement has been previously discussed theoretically [|J , the treatment presented in that 
previous discussion needs to be modified in the case of a complex-valued susceptibility, as 
has been pointed out by Ma et. al. 0. 

The remainder of this paper is organized as follows. Section II describes the extension 
of the spectral theory to uniaxial polycrystalline materials. The application of this theory 
to the nonlinear response of polycrystals is given in Section III, followed by a numerical 
example in Section IV. 

II. SPECTRAL THEORY FOR THE EFFECTIVE MACROSCOPIC LINEAR 

RESPONSE OF A POLYCRYSTAL 

We consider a polycrystalline dielectric material characterized by a position-dependent 
uniaxially symmetric dielectric tensor, which we express in the form 

e(r)=i?- 1 (r)eli?(r), (1) 

where 



?d = e 2 (2) 




is the dielectric tensor in the frame of principal axis and R(r) is a position-dependent or- 
thogonal matrix characterizing the microstructure of a particular specimen (specifically, it 
describes the local orientation of the principal axes with respect to the lab axes). If the 



sample is macroscopically isotropic, it is reasonable to assume that on a large scale its di- 
electric behavior can be characterized by a scalar dielectric constant e e . e e may be defined 
by the relation 

D = -|D(x)rf 3 x = e e Eo. (3) 

where E = y JE(x)<i 3 x is the space-averaged electric field. We assume that E is real 
and is directed along the z axis: E = E z. In general, the fields E and D (as well as the 
dielectric tensor) are represented as complex quantities; the physical fields are related to 
them through E phys (x) = Re (E(x)e- iw4 ), D phys (x) = Re (D(x) e -^ 4 ). 

In the quasi-static approximation, the electric field is given by the negative gradient of a 
scalar potential. We may express this potential through the relation by E = —E V<f), where 
E is real and is the solution to the following boundary-value problem: 

V • (e (x)V0) = in V, 

cj) = cj)Q = —z on S, 

where S is the boundary surrounding V. Using the boundary conditions for <fi and the 
Maxwell's equation V • D = 0, we can show that ||: 

^ J D ■ EdV = i J D • (-E V<f))dV = i J V ■ [D(-£ o 0)] dV 

= 1 fn(-E 4>) ■ dS = 1 y D(-£ O O ) • dS (5) 



i J V • [D(-£ o 0o)] ^ = 7/ / D " {E Vz)dV = Do ■ E . 



Therefore, the definition (Rf) is equivalent to 
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^7 £ n 2 



The result (6) is the equation we use below to express e e in terms of eigenvalues of a linear 
operator. 

In order to achieve this reduction, we first note that instead of the position-dependent 
tensor e* we can use 



- = 1 - uR^CR, (7) 

^2 



where the parameter u is defined by 



«=1-^, (8) 

^2 



and C is a matrix with C*n = 1 and other components equal to zero. We can use this result 
to rewrite the first line of M) as 



V 2 = u (vjR- 1 )) (RV^t^udiRuRydrf, (9) 



where we have used (R 1 J=(i?J, and also employed the convention that repeated indices 
are summed over. From eq. (||), we see that (|j) is equivalent to the integral equation 

(f) = - z + u r<j). (10) 

Here the linear operator T is defined by its effect on a function through 

V(j) = - fd 3 r'G(r,r') (VR- 1 (r')) i (^(r')V'0(r')) i (11) 

and G(r, r') is a Green's function for the Laplace operator: 



V 2 G(r, r') = -<5 3 (r - r')for r in V 
G = for r on the boundary 



(12) 



It is now convenient to define a scalar product of two functions by 

(01-0) = JdV (v<t>*R~ 1 ) i (HVV) 1 , (13) 

With this definition, we can show that T is self-adjoint, non-negative, bounded linear op- 
erator. To show the self-adjoint property, we integrate ( |TT1) by parts using the boundary 
conditions for G to obtain 

T<j> = JdV'(%G(r,r')) R li (r')R lj (T , )d' j <l>(r'). (14) 

Then, using the fact that G satisfies G(r, r') = G(r', r), we find that 

(0|rV> = jdV j dV'd^RuRijdj (d' k G{r, r>)R\ k R\ 1 &^ (r')) 

= JdvJdV'R H R lj R' lk R' ll d j (^G(r,r'))^0(r)^(r') (15) 

= (rv#>. 

To prove that T is real, bounded and non-negative, we consider the eigenvalue problem 

T0i(r) = Si(f>i(r) for r in V, , . 

(pi = for r on the boundary, 

where Si = 1/itj and Ui is the value of u at one of the eigenstates fa (the so-called "electro- 
static resonances"). The physical significance of the latter has been discussed elsewhere [0. 
Next, we note that the problem defined by eq. ([16]) is equivalent to the problem 



V 



(sil - R^CR) Vfa 
at the boundary. 



(17) 



as can be seen by comparing the steps going from (|J) to (|T0D. But from eq. (|H]), we can 
write 

= j dV(j)*V ■ (tul - Rr l CR) Vfa 
= - fdV (^IV^I 2 - \(RVfa) x \ 2 ) , (18) 



from which it follows that < s» < 1. The limiting case Si = 1 (ei = 0) could be realized 
only if the tensor R were position- independent. But this would lead to e e = *e(r), i. e., 
a tensor, which contradicts our assumption that e e is a scalar value. The fact that the 
eigenvalues of (|T6|) are limited to the semiclosed segment [0, 1) proves our statement that T 
is real, bounded, and non-negative. 

From the properties of T, we conclude that the eigenfunctions \4>i) of flT6|) form a complete 
orthogonal set with respect to the scalar product (|Tj|). Hence, the solution to eq. ( PH| ) can 
be expressed (except on the boundary) as 



W = ( ^-i)- Ws5: (_!_)l^!£>. (i 9 ) 

t\Si-s/ {(pi\(pi) 

We can now use this equation to find an analytical representation for e e . 

We begin by using Green's theorem, the boundary conditions in (|j) for 0, and the 
Maxwell equation V • D = to rewrite eq. (|J). The result is 

— = — 777^ I (-E V(f>) ■ BdV = — \— lzT)-dS = — — /z-DdV 
e 2 e 2 VElJ v v> e 2 VE J e 2 VE J 

-1 /".//- « ^ti\ _ a „. 1 /" . ,„ u 

If we now introduce a function 

F{ 8 ) = 1-^, (21) 

^2 



fz-((l- uR- l CR) V(f) dV = — I zz ■ dS + ^ f (VzR- 1 ) i (Sv^) i dV (20) 



then on substituting (|T^) we find 

u 1 1 

F{s) = --{z\<P) = -{z\--^\z) 

i v IWi>| 2 / i \ m) 

Vi (<f>M \s-8iJ' y ] 

This final result is identical in form to Bergman's expression for the analogous function in 
scalar composite materials. The only difference lies in the definition of the scalar product 
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III. APPLICATION TO THIRD ORDER NONLINEAR RESPONSE OF 

POLYCRYSTALS 

As has been suggested by several authors (see, for example, 0), the nonlinear suscep- 
tibilities of composite materials may be hugely enhanced by large fluctuations in the local 
electric field in these materials. The basic idea is as follows: since these nonlinear suscep- 
tibilities depend on higher powers of the local electric field than does the linear dielectric 
function e e , any enhancement of that field will produce an even larger enhancement in those 
susceptibilities than in e e . 



A theoretical description of this enhancement has been given by several authors, initially 
for isotropic composite materials ||, and more recently for polycrystals ||||. The original 
exact expression given in |J is generally difficult to evaluate without approximations. One 
useful approximation involves a decoupling assumption: a certain average of the fourth 
power of the electric field, which enters the exact expression, is approximated as a product 
of averages of second powers ||. However, this decoupling approximation (as well as the 
original exact expression) must be modified slightly when the material of interest has a 
complex-valued dielectric tensor. The need for such a modification was first noted by Ma et 
al, who also generalize the approach of |5| for the case of components with complex scalar 
dielectric functions. In what follows, we further generalize the approach of Ma et al || to 
the case of polycrystals, using the results of Sec. II. 

We consider a polycrystalline material in which D(x, u) and E(x, u>) are related by 

A = CijEj + XijkiEjE k Ei, (23) 

where we suppress the frequency and position dependence of all quantities and sum over 
repeated indices. Next, we assume that a sufficiently large sample of this polycrystal can 
be treated as macroscopically isotropic. Thus, the effective response at the fundamental 
frequency u is given by 

Do = (D) = e e E + x (Eo ■ K) Eo + X (E„ ■ E„) E*, (24) 

where E is the spatial average of the electric field. In component notation, this may be 
written 

Do,i = £ e Eo,i + xfiijfikiEQjE 0tk E j + x^jSmEojEo^Eq^, (25) 

where (...) denotes a volume average and D Qi and E ^ are the i th components of D and 
E . The method of ref. 0] does not permit the two effective susceptibilities to be easily 
calculated independently, but their sum is readily computed. Generalizing eq. (13) of ref. f| 
to the case of finite frequencies, we obtain 

x ^ x + i= (x MM E 1 E i E i ET t (26) 

A) 

where E t = -Ej(x, u) denotes the Cartesian component of the local electric field at frequency 
uj in the corresponding linear polycrystal. 

In this paper we will assume that the fourth-rank tensor Xijui^i u) has certain symmetry 
properties which cause many of its components to vanish. Specifically, we will assume 
that the only non- vanishing components of Xijklfau) (in a frame of reference where the 
coordinate axes are parallel to the local symmetry axes of the crystallite) are those such 
that the indices are equal in pairs. Then Eq. (f26|) takes the form 



Xe — Xiijj ^4 5 \^i) 

where E t = -Ej(x) is the field component parallel to the i th principal axis at x (where we 
suppress the frequency index uj). 
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We approximate the right-hand side using the nonlinear decoupling approximation 
(NDA) fUH, which is specified by the assumption 



;^i^r>«(£?)(i^i 2 ). (28) 

Then using the expression e e = y^ Z)f=i e« / .^(x) 2 ^ we immediately get ||: 

Here the partial derivative denotes 9e e (ei, 62, e^/dei. In the case of a uniaxial material, one 
should calculate this derivative first with e 2 7^ 63, and only then take the limit e 2 = £3. 

The second average on the right-hand side of (^) can be evaluated |§ with the help of 
the spectral approach developed in the previous section. From (|T9|) and our definition (|13D 
of the scalar product, we find 

{lEl] ] ~ T m -v^( ai -^ at -s)(m){m) {<l>M 

£ 2 V \s\ 2 mz)\ 2 ,„ m 

" v Vh-*l 2 (^il^*) ' 

where we have used the orthogonality condition (<fii\<fij) = Sij. 

To evaluate (\E 2 \ 2 ), we note that the boundary conditions (§) for <fi and for 0* are the 
same, since <$$ is real on the boundary. Therefore, we can use E* in place of E in the 
transformations ([5]), so that the definition @ becomes 

e e £ 2 = 1 y E* • DdV = eid^il 2 ) + 2e 2 (\E 2 \ 2 ), (31) 

where we used (I-E2I 2 ) = d-^l 2 )- Hence, 

U ' ' 2 l e 2 e 2 £ 2 / ° 



= i(l-F( a )-(l-l/ S )^ij^ (32) 

if, ^ w - **) m\z)\ 2 \ F 2 

Given an approximation for the effective linear response function e e , the above formulas 
allow us to calculate the enhancement of the non-linear susceptibility. 

The simplest approximation for e e ( e i> e 2,e3) is the effective-medium approximation 
(EMA) [g, which gives 

E^ = 0, (33) 

or for a uniaxial material, 



ei + 2e e e 2 + 2e 
Then the function -F(s) [eq. (21)]is given by 



ei - e, + 2 JiZii. = „ (34) 



F W = 5( 1 _,/iziZ£ | . (35) 



The corresponding T operator has a continuous spectrum, so that the sums in (JZ2"|), ( 50 ) 
and (|32|) should be replaced by integrals. From 

F = /- 1 K^_ dx (36) 

JO s — X 
we find that 



^( z ) = --im[F(x + z0)] = j_ / (x 8 / 9 ) g( _ x) g (x _ g/g^ (37) 

where ^(x) is the usual step function, i. e., 9(x) = 1 for x > and ^(x) = for x < 0. In 
order to evaluate the effective nonlinear response, this expression should be substituted into 
the integrals 



(\E 2 f) = -Jl- dx ^\-^^> El (39) 



^ - 1 -(l-[ 1 dx { ^^f Xy ^- 
2\ Jo |s-x| 2 , 

One immediate consequence of the EMA is that the integral in (^9j) diverges as s — » 
(i. e., as ei/e 2 — > oo), and hence, (I-E2I 2 ) also diverges in the same limit. This divergence is 
related to the divergence of (E%) [eq. (29)]; the physical origin of that latter divergence was 
discussed in Q). Thus, in the NDA/EMA approximation, if X2222 7^ Xe becomes arbitrarily 
large when ei/e2 — > 00, both at zero and at finite frequencies. 



IV. NUMERICAL EXAMPLE 

To illustrate this discussion, we consider a simple model for a polycrystalline quasi- ID 
conductor. In the high-conductivity direction, we assume a Drude metal with dielectric 
function 

eiM = l- , j-, v (40) 

uj{uj + i/r) 

In the perpendicular directions we assume a constant dielectric function 

e 2 = e 3 = 1. (41) 

8 



The resulting complex frequency-dependent e e (u), as given in the EMA, is shown in Fig. 1. 
Fig. 2 shows the corresponding NDA prediction [eqs. (P7|)-(32)] for the enhancement of 



the cubic nonlinearity in the same polycrystal, under the assumption that e e (and F(s)) 
are given by the EMA. For comparison, we show in Fig. 3 the results of an earlier and less 
accurate approximation [|J. The results in Fig. 3 are obtained by using \(Ef)\ instead of 
(\Ei\ 2 ), and evaluating (Ef) using eq. (p9|) . 

Quantitatively, based on our simple example, the effects of the correction noted by |3| 
appear to be relatively minor if only the diagonal elements Xuu 7^ 0, but are more substantial 
if X1122 7^ or X2211 7^ 0. Nonetheless, the corrections do introduce a nonzero correction to 
all the elements of the tensor \. The numerical correction can be seen always to increase 
the absolute value of the corresponding matrix element of x- This trend can be qualitatively 
understood as follows: the spatial average (Ef) which enters into the uncorrected matrix 
elements can, in principle, even vanish under some conditions, but the absolute value (\Ei\ 2 ) 
can never vanish. 



V. DISCUSSION 

Next, we briefly discuss the limitations of the present approach, and the validity of 
the approximations made. At various points in this paper, we have made the following 
approximations: 

• the quasistatic approximation; 

• the nonlinear decoupling approximation; and 

• the effective-medium approximation. 

We now discuss the limitations of each of these approximations. 

The quasistatic approximation is embodied in eq. (£|), which implies that the electric field 
can be expressed as the negative gradient of a scalar potential. This assumption is still valid 
at finite frequencies, provided that the material of interest lies in the long-wavelength limit 
(see, for example, ref. ||). A polycrystalline material is likely to fall in this long- wavelength 
regime, provided that the typical size of a crystallite is small compared to the wavelength 
of electromagnetic radiation in the medium. At optical frequencies, this condition requires 
crystallites of linear dimensions only a few hundred A, but the same approximation might 
hold at microwave frequencies even for micron-size crystallites. More generally, any struc- 
tural correlations in the polycrystal should exist only a scale which is small compared to 
the wavelength; otherwise, there is likely to be significant scattering of electromagnetic ra- 
diation and the quasistatic approximation will break down. Our derivation of the spectral 
representation for a polycrystal is based on the quasistatic approximation. Furthermore, our 
definition of an effective linear dielectric function e e [eq. @] presupposes the quasistatic ap- 
proximation. Indeed, if the quasistatic approximation does not hold, the composite cannot 
easily be described in terms of an effective dielectric function. 

The nonlinear decoupling approximation (NDA) is a way of approximately calculating 
the fourth moment of the electric field in the polycrystal, by breaking this up into a product 
of two second-moment terms. The NDA is known to be quite inaccurate near percolation 

9 



thresholds in conventional composite materials and most likely also in polycrystalline ma- 



terials |10|JTT| . The reason for the inaccuracy is that the NDA neglects local fluctuations in 
electric fields which become very important near a percolation threshold. However, if one 
does not make the NDA, than there is no easy way to express this nonlinear susceptibility 
in terms of the spectral function which describes the linear properties of the polycrystal. 

Eqs. fl27|-J30|) and fl32|) , which are based on our use of the quasistatic approximation and 
NDA, are equally valid for any microgeometry of a polycrystal and can be used with any 
desired approximation for the second moments. On the other hand, the spectral function 
F(s), corresponding to the actual distribution of the electric field, may be sensitive to 
the particular arrangement of the crystallites. We use the effective-medium approximation 
(EMA) to calculate the spectral function which characterizes the linear dielectric function e e . 
Although the EMA does predict the occurrence of a percolation threshold in a polycrystalline 
material, the approximation is likely to be quite inaccurate near that threshold, since it treats 
each crystallite as being embedded in an effective environment. Note that the spectral 
representation itself is more general than the EMA, since it is always applicable in the 
quasistatic approximation. Hence, the spectral representation can be used in conjunction 
with other, more accurate methods of calculating the linear response e e , which take better 
account of local environment of a given crystallite, if such methods can be found. 

We emphasize again that neither the NDA nor the EMA are necessary approximations; if 
better approximations for the fourth moment and for the linear response are available, then 
these can be used to compute the cubic nonlinear susceptibility of a polycrystalline material. 
The inaccuracy of the NDA and the EMA is partially compensated by the simplicity of these 
approximations, which allow many properties to be computed nearly in closed, analytic form. 

To summarize, in this paper, we have extended the spectral representation of Bergman so 
that it applies to the linear effective dielectric function of a uniaxial polycrystal. The exten- 
sion is straightforward, but should be useful in a wide variety of materials. As an illustration, 
we give the spectral function for a polycrystal in which e e is given in the effective-medium 
approximation. Finally, we use this spectral function to calculate the cubic nonlinear sus- 
ceptibility tensor \ e for a uniaxial polycrystal in the nonlinear decoupling approximation 
(NDA), once again calculating the required electric field averages within the EMA. As in 
two-component composites of isotropic materials, the expressions for \ e are slightly altered 
from previous results when one properly accounts for the fact || that the averages (E 2 ) 
and (|-E| 2 ) are unequal. We also give a brief discussion of the conditions under which these 
various expressions and approximations are applicable to real polycrystalline materials. 
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FIG. 1. Re(e e ) (bold line) and 100Re(a e ) (light line), as given by EMA for a polycrystalline 
sample of a quasi- ID conductor. The single-crystal dielectric tensor is assumed to have principal 
values given by Eqs. (pQ) and (||l|) with oj p t = 10. a e is defined by a e 
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FIG. 2. (a) Real and imaginary parts of Xe/Xiiii (bold and light curves) for a polycrystalline 
material, calculated under the assumption that the only nonzero component of the single-crystal 
nonlinear susceptibility tensor is xnii (axis 1 parallel to the high-conductivity axis). The calcula- 
tions are based on EM A results for the linear response (see Fig. 1), and on the model single-crystal 
dielectric tensor assumed in that Figure, (b) Same as (a) except that Xe/xii22 is plotted, assuming 
that only X1122 is nonzero. [The same plot will describe the enhancement of X1212 and X2ii2 ; as 
follows from the definition (p3|).] (c) Same as (a) except that Xe/X22ii is plotted, assuming that 
only X2211 is nonzero. (The same plot will describe the enhancement of X1221 and X2121O (d) Same 
as (a) except that Xe/X2222 is plotted, assuming that only X2222 is nonzero. 
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FIG. 3. Same as Fig. 2, except that instead of the more accurate formulas used in that Figure 
we use the expressions given in Ref. Hj. Plots (a)-(d) correspond to the plots (a)-(d) of Fig. 2. 



14 



